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Abstract 

We present some evidence that in lattice gange theory the condensate contains a 
non-pertnrbative contribntion proportional to A^, the sqnare of the physical scale. 
This resnlt is based on an analysis of the Wilson loop plaqnette expectation from 
Monte Garlo simnlations and its pertnrbative expansion compnted to eight loops. 
The analysis is not fnlly conclnsive since the calcnlations are done on a hnite 
lattice and one needs an extrapolation to inhnite lattice. It has been recently 
snggested that in the glnon condensate a A^-contribntion conld be present coming 
from the large momentnm behavionr of the rnnning conpling and not connected 
to operator prodnct expansion. 


^Research supported in part by MURST, Italy and by EC Programme “Human Capital and Mobility”, 
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1 Introduction 


The Wilson [|^ operator product expansion (OPE) has been applied to the study of 
non-perturbative contributions of physical observables in asymptotically free theories. In 
these studies an important role is played by the gluon condensate (osTr F^). This quantity 
depends on the regularization scheme and on a momentum scale Q, a subtraction point or 
an ultraviolet (UV) cutoff. According to OPE the gluon condensate has the expansion 

tv = + (AVO*) 1^4 + ■ ■ ■ , (1) 

where A is the physical scale of the theory related to the running coupling Oj = as(Q^) at 
the scale Q. For large Q one has Is?/Q'^ ~ exp{—l/47r/3oas}, with the first beta function 
coefficient. Since A has no expansion in as, perturbative contributions are present only in 
Wq, which from power counting is quartically divergent in the UV region. The term with the 
power h?/Q"^ is absent since there are no gauge invariant operators of dimension two. The 
term Wj is the “genuine” vacuum expectation of the gluon condensate of dimension four. 

The fact that the perturbative expansion of W in as is non-convergent |Q|, implies that 
the various terms Wq, IT 4 , ■ ■ ■ cannot be identified simply by collecting contributions with 
a given power of A. Indeed one expects that Wq itself contains non-perturbative terms 
proportional to A^/Q^. These contributions, needed to regularize the non-convergent expan¬ 
sion in tts, have a simple origin. According to the ITEP formulation, Wq is obtained from 
contributions of Feynman diagrams in which the virtual (Euclidean) momenta k are in the 
perturbative region, i.e. ^ A^. The presence of this lower integration bound introduces 
in Wq terms proportional to A^/ Q‘^ which should be canceled by similar terms present in 
the coefficient of the unit operator contribution in the OPE of a physical observable. 

In this paper we address the question whether the non-perturbative contributions in Wq 
are only the ones discussed above which are proportional to A^/Q^. The gluon condensate 
is particularly suitable for discussing this question since no contributions proportional to 
h?/Q"^ are predicted by OPE. Then one may ask whether a non-perturbative term of order 
I?jQ^ could be present. We perform the analysis in lattice gauge theory in which, due to 
the presence of an UV cutoff, the gluon condensate has no UV renormalon [Q 

By using the results on W in SU{3) lattice gauge theory obtained by Monte Carlo 
simulations and by higher order perturbative calculations @ we present some evidence 
that in the lattice regularization W contains terms proportional to Is?/Q‘^. The analysis is 
not fully conclusive since the calculations are done on lattices of finite size and one needs an 
extrapolation to infinite lattice. 

If the presence of terms of order A^/Q^ in W is confirmed, one may worry that OPE 
could be violated. 

Recently Grunberg and Akhoury and Zakharov [g pointed out that the gluon con¬ 
densate could contain a h?/Q^ contribution not accounted for by OPE but originating from 
from power corrections in the high frequency part of the running coupling. Since the per¬ 
turbative Feynman diagrams for the gluon condensate are quartically divergent in the UV 
region, power contributions in the running coupling, which are highly subleading at large 
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momentum, could generate a term in Wq. Power corrections in the running coupling 

at large momentum are naturally expected in physical schemes for the coupling dehnition 


such as the dispersive method 
given in . 


and the V-scheme 0. A similar observation has been 


The paper is organized as follows. In Sect. 2 we recall the expected behaviour of the 
perturbative coefficients for Wq, a quartically divergent quantity. In Sect. 3 we recall lattice 
gauge theory results on the high order perturbative coefficients of W and recall that they 
are consistent with the factorial growth predicted for a quartically divergent quantity. In 
Sect. 4 we analyze the remainder of the expansion and show the indication for the presence 
of a A^/Q^ contribution in the gluon condensate. We show that hnite size effects should not 
be too important. In Sect. 5 we discuss this result and its compatibility with OPE. 


2 Factorial growth of perturbative coefficients 


In the lattice theory all frequencies are bounded by the UV cutoff Q = n/a with a the lattice 
spacing. The condensate W can be written in the general form (we assume an infinite lattice 
for the moment) 


W 


fQ" k‘^dP 
0 


fikVA^), 


( 2 ) 


where k is the softest virtual momentum. This expression is based on the fact that the 
associated observable has dimension four and is renormalization group invariant so that 
the function /(fc^/A^), which does not depend on Q, for large Q, can be expressed in terms 
of a running coupling at the scale k!^. 


We now consider the contribution Wq which is given by the large frequency part of this 
integral. At large the running coupling can be approximated by the perturbative coupling 
as(fc^) obtained by taking into account the hrst few terms of the beta function. For instance 
at two loops one hnds as(fc^) given in terms of the physical scale A by 

/ 1 \bl/bl 

TT- bo = n/{WY, 6i = 102/(47r)^ (3) 

yAnboasik^) J 


For simplicity we may limit the discussion to the case in which at large k'^ the function 
/(/c^/A^) is proportional to the perturbative running coupling. We then consider the contri¬ 
bution to Wq given by 





de 


as{k‘^), 


( 4 ) 


with p S> 1 to make as(fc^) within the perturbative region. Terms with higher power of 
as(fc^) will give similar contributions (see later). depends on the parameter p which 

sets the separation of the low and high frequency part of Eq. (^. Clearly the lower bound 
pA^ contributes to this function by terms proportional to 


Notice that, since the integral is quartically divergent for > oo, the precise form of 
as(fc^) at large k‘^ is very important. In the following we assume that at large the running 
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coupling as(/c^) is given by the beta function at two loops. We introduce the variable 
= 2:0 (l - as/«s(fc^)) , «s = as(Q^), ^ ^ 

and by using for asik"^) the two-loop form one hnds 

— asik'^) ~ dz {zq - z)~^~^ , dTra^ = 6 /^ , 7 = 2 ^, ( 6 ) 

V Oq 

where higher orders in as(fc^) have been ignored. We have introduced the inverse coupling 
fd for the comparison with the lattice theory. The integration region in (^) is mapped into 
the region 0 < z < zo_ = zq (1 — P/P) with Q/P = AnasipA.'^), and one hnds the regularized 
“renormalon” expansion 

PZQ_ 

WJ"" = Af ~ dz (zo - = y p-^ {cr + 0 {e-^y} , ( 7 ) 

Jo 

with e~^°J ~ A^/and cy the renormalon coefficients 

cr = j\f'r{e + yzpr ( 8 ) 

In the coefficients of P~^ in Eq. (J^) the non-perturbative corrections of order hP/QP" are 
essential to make the expansion convergent. They are coming from the lower cutoff ptP 
for the perturbative region and thus they depend on the parameter p. The perturbative 
coefficient grows factorially due to an infrared (IR) renormalon of dimension four. 

A similar factorial growth of the perturbative coefficients is found if one considers the 
contributions from higher powers of the coupling. Then the expression gives a general 
form of the perturbative factorial growth with the numerical constant M which takes into 
account higher order corrections. 


3 Summary of SU{3) lattice gauge theory results 


We hrst recall the basic element of SU{3) lattice gauge theory and then the results on the 
high order perturbative coefficient calculations. The action is 

^[f/] = ^ Tr {Up + f/t), 47r«f = 6/Aat, (9) 

D p 

where the sum extends to all plaquettes P in a hypercubic lattice in four dimensions. The 
plaquette held Up is obtained from the link variable Ufj_{x) = exp{a A^(a:)/A/diat}) given by 
an exponential map on SU{3), where a is the lattice spacing. In the continuum limit a —0, 
at the classical level, the plaquette tends to the Lagrangian density and Eq. (H) tends to 
the Yang-Mills continuous action. In the lattice regularization one has an ultraviolet (UV) 
cutoff Q = vr/a, and is the coupling at this scale. If the lattice is hnite one has also an 
infrared (IR) cutoff Qq = 27r/Ma, with M the number of lattice points in each direction. 
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The condensate is given by the expectation value of the elementary plaquette 

ITixi = 1 - i (Tr Up) . (10) 

We consider also the condensate 1^2x2, the expectation value of the double plaquette, that 
is a Wilson loop around a square of size 2. We denote in general by W these quantities. 

W can be obtained numerically by Monte Carlo simulations. The lattice is finite and 
we denote by W{M) the quantity computed on a hnite lattice of size M. For the hnite 
lattices considered in the present simulations the values of /3iat are taken into the range 
Aat = 6 — 6.5. In this region one expects that the renormalization group properties are 
satisfied: the hnite lattice size is not crucial for /3iat ^ 6.5; the lattice discretization artifacts 
are small for /3iat ^ 6. 


3.1 Factorial growth of lattice perturbative coefficients 


The coefficients cY^{M) of the perturbative expansion of W{M) on a lattice of size M 


/3] 


lat ) 


( 11 ) 


e>i 


are known up to eight loops. The hrst three terms have been computed analytically |]K 
for an inhnite lattice. Eight terms for the expansion of both kFixi and kF 2 x 2 have been 
computed numerically in Ref. [^] for a lattice with M = 8. 


In Ref. P] it has been shown that the growth with i of the hrst eight coefficients is 
consistent with the factorial behaviour described in the previous section for a quantity of 
dimension four like W. In the following we recall the analysis done in p, |T^. 


In order to conhrm that the computed eight coefficients growth factorially as the renor- 
malon coefficients in (H) one has to take into account two facts: i) the numerical calculations 
are done on a hnite lattice; ii) the coupling in (|^) and in are not necessarily within 

the same regularization scheme. 


i) Finite size lattice. The numerical coefficients have been computed on a hnite 

lattice with M = 8 points in each direction. The ehect of the presence of a hnite volume can 
be estimated by putting the IR cutoh Qq = 2ti/M a in Eq. (|), i.e. z < z-^t- = 41n(M/2)//3 in 
(0). This reduces the size of the perturbative coefficients and, for values of [3 with Zir < Zq, 
makes the integral well dehned and the perturbative expansion convergent. We then dehne 

/'min(zir, zo_) 

=M ~ dz {zo - z)-^-'^ = y Ur{M) j3-^ , ( 12 ) 

•^0 ^>1 


where are given by incomplete Gamma functions. For M 

and then c^®“(M) tend to the inhnite volume coefficients of Eq. (| 
coefficients for hnite M been studied in detail in [113 


oo one has Zir —>■ oo 
The reduction of the 
The conclusion is that even for a small 


size lattice with M = 8, the factorial growth up to £ = 8 is not tamed (see also later). 


ii) Continuous and lattice coupling. The two coupling of eqs. oni are not necessarily 
within the same regularization scheme. Diherences in the couplings do not modify the 


4 










factorial growth, but changes subasymptotic contributions. In Q we have assumed that 
they are perturbatively related at two loop by 


^ = Aat - r- - 


A. 


(13) 


at 


It is known that the relation between the lattice and continuum coupling involves large 
perturbative corrections. For instance if (3 is the MS coupling one hnds [|^ the values 
r = 1.8545 and r' = 1.667. By using (^) we can connect the coefficients in (|^) with the 
ones in (0) by computing the new coefficients CY'^{r,r', M) 


wr(M) = criM) n 

£>i 




= i: Cl 

i>i 


(r, r', M) Aat • 


(14) 


We find that for large ^ the numerical coefficients of the single and double plaquette agree 
with , M) for large values r = 3.1 and r' = 2.0. See Fig. 1. The fact that the 

resulting value of r is larger than the one coming from MS can be simply explained by 
assuming that the scale entering in the running coupling in Eq. (^) is sk‘^ < in analogy 
to the exact result of the sigma-model (see also |]I^ ). 
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Figure 1: Coefficients of the perturbative expansion as function of the (loop) order The 
open squares represent A®"* from Ref.6 for i < 8. The open tringles represent the renormalon 
coefficients M) in Eq. (14) for r = 3.1, r' = 2.0 and M = 8. 

Therefore the main conclusion of Ref. was then that the factorial growth of the com¬ 
puted coefficients A®"* is in agreement with a renormalon associated to dimension four, i.e. the 
the gluon condensate dimension. 
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4 Evidence of a /Q^ contribution 


Here we present the evidence that in lattice gange theory the condensate contains terms 
of order ~ g-^o/3/2 analyze W — Wq as a fnnction of p. The contribntion W is 

obtained by the Monte Carlo simnlation. The contribntion Wq is constrncted by adding to the 
computed eight-loop perturbative terms a remainder. Since for large orders the perturbative 
coefficients of approach the lattice ones (see Fig. 1), we use this function to obtain 

the remainder. In Eq. (JT^) we assumed that the running coupling as{k‘^) has the two- 
loop asymptotic behaviour. Under this assumption the constructed Wq contains, beside the 
perturbative terms, power terms of order ~ In this case the difference W — Wq 

should be of order A^/Q"^ ~ We shall hud instead that there is a contribution of order 

AVQ2 _ g-^o/3/2_ 

First we plot in Fig. 2 the quantity 

Aiir(M) = W(M) - (15) 

1=1 

in the range j3\s± = 6 — 7 for various values of L < 8. The quantity W{M) is obtained 
from the Monte Carlo simulation of Ref. on a lattice with M = 8 as for the perturbative 
coefficients. 

We observe that in this range of /diat the quantity AlW{M) approaches for L —>• 8 the 
behaviour of A^/Q^ instead than the expected behaviour of 

Before drawing a dehnite conclusion on this result one needs to analyze whether the h? / 
behaviour is modihed by considering: 1) the remainder of the perturbative expansion; 2) the 
effects of hnite volume. 


1) Remainder of the perturbative expansion. As previously recalled, at large ^ the perturba¬ 
tive coefficients r', M) in Eq. (JT^) reproduce the lattice coefficients We then 

estimate the remainder SWq by subtracting from hFo®’^(M), given in Eq. (|^), the hrst eight 
terms in (pAj) 


W„(M) = CT(r,r',M)0,-t- 

We than obtain the following estimate for lUo 


( 16 ) 


£=1 


W„(M) = Y: /3,-J + m(M). ( 17 ) 

e=i 

We plot in Fig. 2 the quantity 

AhF(M) = W(M) - Wo(M), (18) 

for M = 8 in the region /diat = 6 — 7. We see that the behaviour hf is still maintained. 
The conclusion is that in the region considered for j3\st the hrst eight perturbative terms give 
a reliable approximation of Wq, at least for M = 8. 


2) Finite volume. This effect is quite difficult to estimate without performing a direct Monte 
Carlo simulation on lattices with M sufficient large to have the IR cutoff below the Landau 
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Fig.2a 


Fig.2b 
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Figure 2: (a) The subtracted MonteCarlo data A^W of Eq. (15) compared to and 

for various values of L: upper curve for L = 2, lower curve for L = 8; (b) The sub¬ 
tracted MonteCarlo data AW of Eq. (18) after resummation of the renormalon contribution 
compared to A^/Q^ and 


singularity, i.e. ln(M/2) > (3/12bo. For the contribution Wq we can estimate the effect of 
the hnite volume. For the hrst eight coefficients this has been done in and, as already 


recalled, for M = 8 the factorial growth is still present. We can study the M dependence of 
the remainder 5Wq{M) in (pTf ) and we hnd that in the considered region of /diat the effect of 
hnite size is small, i.e. less than 5% . 


5 Discussion and conclusion 


One has to consider the following two indications. 

1) As shown in [^, the hrst eight perturbative coefficients of W seem to agree with the 
factorial growth corresponding to a IR renormalon associated to an operator of dimension 
four, as required by OPE (see hg. 1). 

2) Here we have studied the contribution Wq in (|^) obtained from the hrst eight terms 
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of the perturbative expansion and a remainder (^) constructed on the hypothesis that 
only corrections are present. By subtracting from W the term Wq we have found 

indications of an additional contribution proportional to (see £g. 2). 


Some caveat are in order. The reason for the unexpected behaviour could be 

that our analysis is not complete. The major problem is the hniteness of the lattice size. 
However we have estimated its effects on the considered contributions to Wq and they seem 
to be small. It may be that the Monte Carlo simulation for W contains spurious hnite 
size effects giving an effective A^/Q^ behaviour. Excluding this possibility would require an 
investigation on a very large lattice with ln(M/2) ^ f3. 


Recently it has been argued by Grunberg and by Akhoury and Zakharov p] that terms 
of order A^/Q^ can be present in the gluon condensate which are not accounted for by OPE, 
but are due to power corrections in the running coupling at high momentum. In physical 
schemes [P, highly subleading power corrections at large momentum are naturally present 
in the running coupling. A similar observation has been given in [^]. These corrections could 
be responsible for the appearance of A^/Q^ terms in the condensate due to the fact that the 
integral for W is quartically divergent. A A^/ k"^ contribution in as{k‘^) in the integral gives 
two terms. The first of order A^/Q^ is coming from the UV region the second, 

of the canonical order A^/Q'^, is coming from the IR region pA^). These A^/Q^ terms 

are of “perturbative” nature and are then naturally associated to the contribution Wq in the 
OPE. Moreover they should be process independent as the running coupling. An important 
question is whether these A^/Q^ are phenomenologically relevant (see |]^, ^). 
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